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Abstract 

The technique of conformal mappings is applied to enlarge the convergence do- 
main of the Borel series and to accelerate the convergence of Borel-summed Green 
functions in perturbative QCD. We use the optimal mapping, which takes into ac- 
count the location of all the singularities of the Borel transform as well as the present 
knowledge about its behaviour near the first branch points. The determination of 
a s (m^.) from the hadronic decay rate of the r-lepton is discussed as an illustration 
of the method. 



1 Introduction 



The behaviour of the large order terms in perturbative quantum field theory has been a 
subject of permanent interest 0-0. Recently, this problem received much attention in 
the case of QCD 110|-|20|. As is known, the creation of instanton-antiinstanton pairs and 



certain classes of Feynman diagrams are responsible for a factorial increase of the large 
order coefficients of the perturbative expansion of the QCD Green functions, making this 
series divergent in the mathematical sense. Moreover, the growth of the large-order coef- 
ficients is so dramatic that, when combined with other difficult circumstances (like their 
non-alternating sign and the extraordinarily small analyticity domain of the Green func- 
tions in the coupling variable ||), it leads to the situation that some of the usually efficient 
summation techniques |2l| are not applicable. One of them, the Borel summation, was 
very much investigated in the recent time. The growth of the large order perturbative co- 
efficients of the QCD Green functions leads, under certain conditions, to Borel transforms 
with singularities in the Borel plane that make the integral defining the Green function 



by the Laplace transform ill-defined. Of course, as discussed in | 19fl , the Borel technique 



is not the only mathematical method by which a divergent series can be summed, but 
an ambiguity emerges in every summation method, once it is discovered in one of them. 
The real source of nonuniqueness consists in a missing piece of information about the 
quantities to be calculated, which adopts different forms in different summation methods, 
but has to be added to eliminate the ambiguity. 

In QCD the Borel nonsummability originates from the infrared regions of the Feyn- 
man diagrams, where nonperturbative effects play also an important role. Therefore it is 
natural to assume that the ambiguities of perturbation theory must be compensated by 
nonperturbative contributions. In fact, it turns out that general concepts like analyticity, 
renormalization group or specific properties of the QCD vacuum are unavoidable when 
discussing the large order behaviour of perturbation theory j|] , || , |19[] . As an example, we 
recall that the argument given by 't Hooft |J for the Borel nonsummability of QCD relies 
on nonperturbative properties, mainly the momentum plane analyticity combined with 
renormalization group invariance. So the perturbative and the nonperturbative regimes 
of the theory cannot be separated, and their interplay is very clear when attempting to 
perform the summation of the large orders of perturbation theory. The Borel plane is 
particularly suitable for discussing this aspect, since the singularities of the Borel trans- 
form offer a very intuitive measure of the ambiguities of the perturbation theory and 
suggest the way to compensate them. The properties of these singularities in some ap- 
proximations (like massless QCD in the large /?o limit, when they are poles [[13]) were used 
recently to provide estimates of the truncation error in the theoretical determinations of 
some accurately measured quantities. 

A natural question is whether it is possible to improve the accuracy of the Borel 
summation using the first Taylor coefficients of the Borel transform known from the 
explicit low order calculations, supplemented with some (approximate) information about 
its singularities in the Borel plane. We address this question in the present paper. We use 
as input the assumption that, for suitable QCD amplitudes, there are no other singularities 
in the Borel plane than those located on the real axis, at a nonvanishing distance from 
the origin. The precise nature and strength of these singularities is not known in general, 
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except for the nearest ones, which can be characterized (at least approximately) by using 
general principles. As discussed in ||, the singularities of the Borel transform require 
the introduction of higher dimensional operators, which ensure the compensation of the 
ambiguities present in the usual perturbative terms by the ambiguities inherent in their 
Wilson coefficients. This allows one also to infer a universal behaviour of the Borel 
transform near the first ultraviolet renormalon | 20fl . On the other hand, as discussed in 
fT0|| , the location and nature of the first infrared renormalon can be plausibly predicted 
too, by nonperturbative arguments. In our approach, we take as input this knowledge 
about the first ultraviolet and infrared renormalons. 

Our purpose was to exploit in an optimal way this information, in order to improve 
the accuracy of the Borel summation in the frame of a specific prescription of handling 
the singularities of the Borel transform. To this end we use the analytic continuation of 
the Borel transform outside the circle of convergence of the Taylor expansion, achieved 
by the technique of conformal mapping. As is known, the conformal mappings are very 
suitable for accelerating the convergence of power series. The existence of an optimal 
expansion, with the largest convergence domain and the best asymptotic convergence 
rate, was proven in 2^] a long time ago. The method is applicable if the position of the 
singularities of the function to be approximated is known or can be reasonably guessed, 
which is the case in many situations in particle physics. 

In the context of the Borel summation in quantum field theory, the conformal mappings 
were first considered in ||-||. More recently, the method was applied in |p3[|-|p4| to the 
Borel transform of QCD Green functions, following a suggestion made in [|Kj. The purpose 
was to estimate (and possibly reduce) the influence of the first ultraviolet renormalon and 
of the associated power corrections on observable quantities, by using a variable in which 
this singularity is pushed further away from the integration range of the Borel transform. 



However, from the point of view of the convergence rate the mapping used in [[23|j -[f24 
is not optimal, since it takes into account only a part of the singularities of the Borel 
transform, the ultraviolet (UV) renormalons. By using an optimal treatment, which 
takes into account also the infrared (IR) renormalons and the behaviour near the first 
singularities, an increased convergence rate and consequently a smaller truncation error 
are to be expected. 

The objective of our work is to establish whether the optimal mapping technique is 
numerically relevant in the Borel plane for situations of physical interest (we mention 
alternative attempts to enlarge the convergence domain of the Borel transform, based on 
Pade approximants [p5] ). To illustrate our discussion we consider, as in |p3fl , the Adler 
function of the massless QCD vacuum polarization and the determination of the strong 
coupling constant a s (m^) from the hadronic r decay rate. In the next section we briefly 
review some properties of the Adler function and of its Borel transform. In section 3 we 
present the technique of optimal conformal mapping and investigate its efficiency in the 
Borel plane on several mathematical models which simulate the physical situation. We 
discuss also the determination of the strong coupling constant a s (m^.) using the present 
technique. Some conclusions are formulated in section 4. 
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2 The Adler function and its Borel transform 



We consider the correlator 

i J d 4 x e iq - x (0| T {V»(x), ^"(0)*} |0) = (g"g" - g^q 2 ) U(s) , (1) 

where s = q 2 and V* = q^q is the current of a massless quark. From the general principles 
of causality and unitarity it follows that the amplitude II (s) is an analytic function of real 
type in the complex plane s, cut along the real positive axis from the threshold 4m 2 of 
hadron production to infinity. It is convenient to define the Adler function 

D(s) = ~Tl(8), (2) 

which is ultraviolet finite and is also analytic in the complex s-plane cut above the unitarity 
threshold. This function was much investigated lately in connection with the determi- 
nation of the strong coupling constant a s (m 2 ) from the hadronic decay of the r lepton 
flPTl, p6 | - [ 29|j . The hadronic decay width, normalized to the leptonic one, is defined 



Rt = r(T 7 V * + hadTOns) = fds , (3) 



as 



T(r — ► u T e u e ) J ds 

where the inclusive hadronic spectrum dR T (s)/ds is related to the spectral part of the 
correlator ([[]): 

dRJs) 3(1 + 5 EW ) / s \ 2 / 2s 



, .1- — 1 + — Imn( S + ze). (4) 
ds irm^. V m^J \ m^J 

The factor Sew — 0.0194 accounts for electroweak radiative corrections. The decay rate 



was measured recently with great accuracy [RT 



Using the analyticity properties of the function IT(s) in the momentum plane, the 
relation (^) can be transformed by a Cauchy relation into 

3(l + fcw) ,d.A_ jW 1 + 4 
2m J s \ mi J \ mi 



Jlr = aiZ2wi^( 1 _. 1+ . )DWi (5) 



where the integration runs along a closed contour in the complex plane, taken usually to 
be the circle \s\ — m 2 . 

The relation @ is the starting point for the computation of the r hadronic width 
in perturbative QCD. At complex values of s the Adler function admits the formal 
renormalization-group-improved expansion 

D(s) = l + f:D n (^l) n . (6) 

„=i V vr 7 

The strong coupling a s (/i 2 ) satisfies the renormalizat ion-group equation 

/i 2 ^^ = -«,(/i 2 ) E PnM» 2 )T +1 , (7) 
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with the first coefficients f3 n defined in terms of the number n/ of quark flavours as 

33 - 2n f 



0o 
Pi 



12vT 

153 - 19n f 
24^ 



The coefficients D n in the expansion (y) were computed for n < 3 p3[- ||37|| . In the MS 
scheme with rif = 3 they are 



^3 



1.63982, 
6.37101. 



(9) 



On the other hand, the large-order coefficients, D n for large n, have the generic factorial 
behaviour 



D„ ?zC k n\n 6 



k 



+ .... 



(10) 



where the index k takes the values: —1, ±2, ±3 



In the Borel method of summation one defines the Borel transform of the Adler func- 
tion as 

oo 

B(u) = y £b n u n , (11) 



n=0 



where 



b n = r 



1 D 



n+l 



Then D(s) can be expressed formally in terms of B(u) by the Laplace transform 



D(s) = 1 + 



duB{u) exp 



u 



(12) 



(13) 



To illustrate our technique we shall use also a Borel-summed expression for the hadronic 
decay rate R T . Such an expression was obtained in |23| by inserting in (|5|) the Borel 
representation (|i~3"D of the Adler function and performing the integration along the circle 
\s\ = ml in the one- loop approximation 



1 



at. -s 



A,ln(- S /A 2 ) 
of the running coupling. This procedure gives |23[| 

\ r<x> 



(14) 



R T = 3(1 + S_ 



EWj 



1 + 



■kPq Jo 



dwexp 



u 



P a s (m 2 T ) 



B{u) F{u) 



where 



F(u) 



-12 sinfTr-u) 



tiu(u — 1)(m — 3)(w — 4) 



(15) 



(16) 
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We shall use dig) as a starting point for a determination of a s (m^) in Section 3. 
The growth ([II]) of the Taylor coefficients Z) n leads to the dominant behaviour 



B{u)^C k T{5 k + l){l 



u, 
k' 



+ ..., 



(17) 



which shows that the function B{u) becomes singular at the points u — k, with k = 
— 1, ±2, ±3.... The precise values of C k and 5^ are not known in general. However, from 
general arguments it was shown that the nature of the first branch points of the Borel 
transform is universal 



T0[| , [p0 |. More precisely, near the first UV renormalon at u 



the Borel transform behaves as 



where 20 



B(u) 



7i 



u 



-71 



Pi 

P 2 o 



3-9 + Al 



(19) 



Here Ai is a parameter depending on the number of flavours, which reflects the mixing of 
higher dimensional operators in the renormalization group equations |2(| . Similarly, near 
the first IR renormalon at u = 2 the behaviour is 



B(u) 



r 2 (2 - u) 



-72 



where 10 



72 



1 + 2^ 
+ V 



(20) 



(211 



Using the first coefficients /3j from @ and the parameter Ai given in |2(J (equal to 0.379 
for rif = 3 and 0.630 for = 5) we obtain 



for rif = 3, and 



71 = 2.589 , 72 = 2.580 



71 = 2.972, 72 = 2.316 



(22) 



(23) 



for rif = 5. We emphasize that only the nature of the first renormalons is known, and 
nothing can be said about the residues T\ and r 2 appearing in (0) and [H], respectively. 

Strictly speaking, the integrals (|T3|) or flT5|) have nothing to do with the summation of 
the perturbative series, because one condition of the Borel theorem (the existence of the 
analytic continuation in the a s plane from the convergence disk to an infinite strip around 
the positive real semiaxis) is not satisfied |2TJ, 0. This can be seen from the singularities 
of the Borel transform given in (|17D: the poles situated on the real positive axis (IR 
renormalons) make the integrals (|l~3|) and (|lj|) ambiguous. In order to compute them a 
prescription has to be adopted, by suitably choosing the integration contour in order to 
avoid the singularities. But this is not the Borel summation. Different prescriptions give 
different results, and a measure of the intrinsic ambiguity of the perturbation expansion is 
given by the difference between these results, if no a priori arguments in favor of a certain 
choice exist. 
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A prescription adopted by several authors ||, (T^|, []17| is the "principal value" (PV), 
defined as 



PV / duexpf-— ) f(u) ' 



2 



/ dwexp ( ) f(u) du + / du exp ( — — ) f(u) 

Jc+ \ aj Jc- V a s J 



(24) 

where C± are two lines parallel to the real axis, slightly above and below it. This definition 
is a generalization to arbitrary singularities of the Cauchy principal value of simple poles, 
and has the advantage of yielding a real result when a s is real. Although this prescription 
is not unique, we adopted it as a working hypothesis. 

We applied the definition (^ip for computing the Borel-summed Adler function (given 
by the Laplace integral (|i~3D), and the hadronic r decay rate (given by (|i~5D). In the first 
case the parameter a s is related to the running coupling a s (— s) (and may be complex if s 
is complex or in the Minkowskian region s > 0), while in the second case it is proportional 
to a s (m 2 T ). We use as input the first Taylor coefficients of B(u) (known from fll2|) and the 
calculated values (^)), supplemented by the knowledge on the location and the nature of 
singularities of this function. As discussed in the Introduction, our purpose is to improve 
the accuracy of the calculation by the technique of conformal mappings, which exploits 
in an optimal way this input information. 



3 Optimal conformal mapping of the Borel plane 

3.1 Remarks on the theory 

The use of conformal mappings for improving the convergence of power series in particle 
physics was discussed for the first time in Ref. |22j| . The problem formulated in was 
to find the optimal conformal transformation which minimizes the asymptotic truncation 
error of a power series, taking into account the location of the singularities of the function 



to be approximated. First we briefly describe the results obtained in [22]. We consider a 
function f(u) analytic in a domain T> of the complex u plane containing the origin, and 
write its Taylor series truncated at a finite order N as 



N 

f {N \u) = Y.Uu n . (25) 



ra=0 



According to general theory, the series (|25|) converges, for iV — > oo, inside the circle 



passing through the nearest singularity of the function f{u) in the complex plane, the rate 
of convergence at a point u situated inside the circle being that of the geometrical series 
in powers of where r = \u\ and R is the radius of the convergence circle. Therefore, 
the convergence rate is strongly influenced by the distance of the singularities of f{u) 
from the origin, and can be improved by using a suitable change of variable, in which the 
singularities are pushed further away from the region of interest. 

Let us consider a conformal mapping w = w(u) of the plane u onto the plane w, such 
that u>(0) = 0, and write the truncated Taylor expansion of the function f{u) in the 
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variable w: 

N 

fL N) (u) = J2c n w n . (26) 

n=0 



As pointed out in |22J, the best asymptotic rate of convergence of this series in a certain 
region of the complex plane is achieved when w is such a transformation of the u plane 
that the corresponding ratio ^ is minimal for every point in that region. According to 



the theorem proven in [22] this is realized by a conformal mapping w(u) which maps the 
whole analyticity domain T> of the function f(u) in the u plane into the interior of a circle 
in the plane w. The proof of the theorem is based on the fact that the circle is the domain 
of convergence for power series, and on the Schwartz lemma, which implies that the larger 
the domain mapped inside the circle, the better is the asymptotic rate of convergence (for 
details see f22|). 



In what follows we shall apply this technique to the Borel transform B(u) of the Adler 
function. The nearest singularities of the function B(u) are situated at u = —1 and u = 2 
and the power expansion (JTTJ) converges only inside the circle \u\ < 1 passing through the 
first UV renormalon. It is easy to see that the optimal conformal mapping in the sense 
explained above is given in our case by 



w= \ S l = — ■ (27) 



By this mapping, the complex u plane cut along the real axis for u > 2 and u < — 1 is 
mapped onto the interior of the circle \w\ < 1 in the complex tu-plane, the origin u = 
of the u plane becoming the origin w = of the w plane, and the upper (lower) lips of 
the cuts are mapped onto the upper (lower) semicircle in the plane w. Particularly, all 
the singularities of the Borel transform, the UV and IR renormalons, are now situated on 
the boundary of the unit disc in the w plane, all at equal distance from the origin. The 
Taylor expansion of the Borel transform in powers of w, 

N 

BL N \u) = J2cnw\ (28) 

n=0 

will converge for iV — > oo up to points close to the renormalons. This is a considerable 
improvement with respect to the usual expansion (|IID , whose convergence in the Borel 
plane is limited by the circle reaching the first UV renormalon. The requirement of 
convergence of (^) on this disc implies the holomorphy of the expanded function on 
the disc. In this way, the expansion in powers of w(u) makes full use of the analyticity 
property that is universally (but tacitly) assumed in all QCD considerations, namely that 
there are no singularities in the Borel plane other than those situated on the real axis, 
at a nonvanishing distance from the origin. This essential, additional assumption has, to 
our knowledge, not been explicitly used. 

For comparison we give the conformal mapping used in p3f- |[24||. 



z= ^~\ (29) 
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which maps the u plane cut along u < — 1 onto the interior of the unit circle in the z 
plane. In the z plane the UV renormalons are situated along the boundary of the unit 



circle \z\ = 1, but the IR renormalons are situated inside this circle. As noticed in f24}| , 
pushing away the ultraviolet renormalons by ( |2"D| ) has a price in moving the first infrared 
renormalon (and actually, the whole positive real semiaxis) closer to the origin. This is 
why the convergence domain of the power series in z is limited by the first IR renormalon 
and, as a consequence, the convergence rate of the series in powers z will be worse than 
that obtained with the optimal variable (p7|) . The use of the optimal conformal mapping 
(P7|) is therefore highly desirable, because it does not suffer from this shortcoming, placing 
all the renormalons onto the circumference of the unit disk. 



As was noticed in [3S|, a further improvement of the convergence rate can be reached 
if some information about the nature of the singularities of the expanded function f(u) 
is available. The idea is that the power variable w(u), taken as a function of u, should 



resemble f(u) as much as our knowledge of f(u) allows it. (As it was put in |22[, if we 
were to know / exactly, the most rapidly convergent expansion would be that in powers 
of / itself, in which case it would reduce to the identity / = /.) 

In practice, however, our knowledge of the expanded function is only approximative. 
For instance, as discussed above, we know that near the branch points u — — 1 and u = 2 
the function B(u) behaves like (u + 1)~ 71 and (2 — u)~ 12 respectively, with the 7, real 
positive numbers. In this case it is convenient to expand the product (m+1) 71 (2— u) 12 B(u) 
in powers of the optimal variable w defined in (p7|) , and introduce then explicitly the 
singular factors. The expansion of the function B(u) will have the form 

1 N 

B ( W N) (u) = r—, — Y c n w n . (30) 

w V ' (u+ l)7i(2 -u)^ ^ n y ! 

The singularities themselves may survive as positive powers (u + l) 71 and (2 — u) 72 , the 
bonus nevertheless being that the positive exponents 7, keep the values of the function 

(u+ 1) 7i (2-m) 72 j B(m) (31) 

finite near u = u\ and u = 112, which softens their numerical impact. This step will imply 
no large-order improvement of the convergence rate (because the rate is given by the 
position of the nearest singularities), unless some of the two singularities is fully removed 
by it. But it may represent a considerable improvement at low orders, even if the nature 
of the nearest singularities is known only approximately. A nice example of efficiency of 
this approach in practice was presented by Soper and Surguladze in [Q. 

3.2 Discussion of mathematical models 



We tested the practical efficiency of the conformal mapping (|27|) for a number of functions 
having logarithmic or power branch points at u — — 1 and u — 2. We took functions close 
to the physical situation as described in Section 2, i.e. we started from a "perturbative" 
expansion of the form 



N 

I 

n=0 



B^ N \u) = Y / b n u n , (32) 
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with low values of N. The expansion (|28| ) in terms of the variable w is obtained by 
replacing u in (|3"2"D with the expansion 



u 



(AT) 



(33) 



n=l 



which follows from the inverse of (|27|), and keeping only terms up to the order N (i.e. 



iu 



(N)\N 



CiW , etc). The numerical values of the first coefficients C n are: 



8 16 40 



224 



81 ' ~ J 243 "' ^ 
For comparison, the expansion in powers of the variable z given in (|29|) is obtained using 



u 



(TV) 



TV 



n=l 



with the numerical values [p3f 

Ci = 4 , C 2 = 8 , C 3 = 12 , = 16 , C 5 = 20 ... 



(35) 



(36) 



We first computed the model functions and their various approximants at points u 
inside the analyticity region, near the origin of the Borel plane. In most of the cases 
investigated the expansions in powers of the optimal variable w approximated the exact 
functions much better than the standard expansion (^) or the series in powers of the z 
variable. This feature was visible even with a few terms in the expansion. Actually, as 
we discussed in the previous Section, for the physical applications we are interested in 
the calculation of Laplace integrals like fll3|) or (|i~5"|) . We evaluated this integral, with the 
generalized principal value prescription defined in ( p4"| ) for a large number of functions of 
physical interest. We consider as an example the function 



B(u) 



+ 



r-2 



(1 + u)ti (2 -it) 



12 



Nib _ 
+ v — 



(37) 



which simulates the contribution of a few renormalons. The principal value (|24j) of the 
Laplace integral was computed numerically with great accuracy. In order to check the 
computations we used the relation f39| 



lim 

e-*o Jo 



00 d ^ exp(-w/a s ) 



(it + b + ieY 



^-" e W^r(-v + l,b/a a ) , Rea s > 0, 



(3? 



where e > and T(u,z) is the incomplete gamma function [39|, analytically continued 
from the region Re z > to the whole complex plane z cut along the negative real axis. 
For integer v this can be expressed equivalently as |40f 



lim 



00 exp(-w/a s ) 



(u + b + ie) r 



-E n (b/a s ) , Rea s > 



(39) 
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in terms of the exponential integral functions E n [z). Actually, as seen from (|37|), in the 
physical case the denominators must be defined so as to give the correct cut structure of 
the Borel transform. This case is obtained from (|39[) as 



e^oJo (\b\-u-ie) 1 



e b/a s 

—E n (b/a s 



~\b\ 



(40) 



where b = —\b\. 

As a side remark, we mention that the above relations are useful for defining the prin- 
cipal value prescription for arbitrary values of a s . First, by means of repeated integration 



by parts in (|39|) we can express the left hand side as |]39 

, \1— ra+m I 



-u/a s ) 



00 exp( 

du- , 

o (u + b + ie) 



n-1 

£ 

m=l 



1)! 



\l-n 



1)! 



+ 



[n 



1)! 



-u/a s ) 



00 , exp 

dM 7 , , 

o (u + 6 + le) 



Rea s > 0. (41) 



We now apply the definition (|24j ) of the principal value and use the symbolic relation 

1 „. 1 



PV- 



=F in5{u + b) 



(u + b±ie) ' ' (u + b) 
in the last term in fl4"I|). We obtain thus the following expression of the principal value 



(42) 



pv fdti^-y 

Jo (u + b) n 



—E n (b/a s ) + m- 



'-as) 



l-n 



1 ! 



Re a, > . 



(43) 



For real values of a s the last term in the above relation is purely imaginary and com- 
pensates the imaginary part of the first term. In this case the definition fl43|) amounts 
therefore to taking the real part of the right hand side of (|39D . For complex a s , when the 
last term in fl43|) has also a nonvanishing real part, the compensation does not occur, and 



the result is complex. As discussed below Eq. (|24j) , complex values of a s appear in the 
Borel summation of the Green functions in the complex momentum plane or in the time- 
like region. For some Minkowskian quantities a definition of the principal value, based 

is general and 



The above expression 



on physical arguments, was proposed in 
covers all these cases. 

In Fig.l we give for illustration the results of our analysis for the model function ( |3"T| ) 
with the parameters: r\ = 1, r 2 = 4, r n = 0,n > 3, 71 = 2.5,72 = 2.5, and a s in the 
range (0.1-0.9). The Laplace integral of the exact function (^7p is indicated together with 
the results given by its ordinary perturbative expansion (truncated at N — 3), and the 
expansions accelerated by the conformal mappings z and w, both in the simple versions 
and with the improvement explained in fl30|). To simulate the physical situation in a more 
realistic way, we assumed that the nature of the first singularities is not exactly known and 
used in the improved version (|3"0D the product of B(u) with the factors (u + 1) 2 '°(2 — u) 



\2.0 



which do not compensate exactly the singularities of the model function (37). 

As seen from Fig.l, the combined technique of conformal mapping and the explicit 
treatment of the branch points, supposing that some (approximate) information about 
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Figure 1: Laplace integral for the model fl3"T|) as a function of a s : exact values (a); 
perturbation expansion (b); expansion in powers of z (c); improved expansion in powers 
of z (d); expansion in powers of w (e); improved expansion in powers of w (f). The series 
are truncated at iV = 3. 



the behaviour of the function near the first singularities is available, improves the accu- 
racy of the Borel integral, especially for large a s . The values a s ~ 0.2 — 0.3 (of interest 
in the hadronic r decay) are on the boundary of the region for which the improvement 
is significant at this order, N = 3. We notice that a major part of the improvement is 
brought by the separate treatment of the branch points, according to (|30|), especially at 
low N. Even the standard expansion in powers of the Borel variable u gives good results 
if the nature of the lowest singularities, assumed to be exactly known, is treated explicitly 
as above. However, when the behaviour near the first singularities is known only approx- 
imately, the expansion in the optimal variable gives in general the best approximation, 
especially when the order N of the truncation is increased. We illustrate this fact in 
Table 1, where we indicate the Laplace integral for a s = 0.3, of the function ( |37|) with 
7i — 72 — 2.5 as a function of the truncation order N, for different types of expansions. 
The "improved" expansions were obtained now by expanding in powers the product of 
B(u) with the factors (u+ 1) L5 (2 — w) 1 - 5 , close but not identical with the actual behaviour 
of (0). 

For larger values of a s the improved accuracy obtained by using the optimal mapping 
is even more impressive. Some results are presented in Table 2 and in Table 3, for a s = 0.5 
and a s = 0.8, respectively. Similar results were obtained also for model functions with 
more singularities on the real axis. 

A closer look at the Tables 1, 2, and 3 reveals that there are essentially three circum- 
stances affecting the convergence properties: (i) the use of a convenient (including the 
optimal) conformal mapping, (ii) explicit (but, in practice, approximative) account of the 
branch point singularities, and (iii) exponential damping of the integrand by exp (—u/a s ). 

The effect of the factor (i) can be seen from the fact that, in each of the Tables, the 
column (c) possesses better convergence properties than the column (a), and the column 
(e) has better properties than the column (c). The effect of (ii) is seen from the fact that, 
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N 


(a) 


(b) 


(c) 


(d) 


(e) 


(/) 


3 


0.35 


0.530 


0.50 


0.546 


0.47 


0.5596 


4 


0.90 


0.610 


0.540 


0.5625 


0.64 


0.5701 


5 


-0.10 


0.5732 


0.547 


0.5722 


0.5457 


0.5701 


6 


2.2 


0.516 


0.5631 


0.5783 


0.613 


0.56518 


7 


-3.6 


0.82 


0.5743 


0.582 


0.63 


0.56518 


8 


13 


-0.12 


0.587 


0.584 


0.515 


0.563102 


9 


-40 


2.3 


0.599 


0.584 


0.64 


0.563100 


10 


143 


-3.4 


0.613 


0.583 


0.503 


0.563470 


11 


-541 


9.5 


0.63 


0.581 


0.540 


0.563467 


12 


2 10 3 


-19 


0.64 


0.5776 


0.590 


0.563783 


15 


-2 10 5 


199 


0.72 


0.56400 


0.5783 


0.563713 


20 


2 10 9 


-8 10 3 


1.1 


0.530 


0.582 


0.563689 


25 


-4 10 13 


3 10 5 


2.7 


0.48 


0.5687 


0.563681 


30 


2 10 18 


-5 10 7 


13 


0.43 


0.559754 


0.563682 



Table 1: Approximations of the Laplace integral for a s = 0.3 for different truncation 
orders: (a) expansion in powers of u; (b) "improved" expansion in powers of u; (c) expan- 
sion in powers of z; (d) "improved" expansion in powers of z; (e) expansion in powers of 
w; (/) "improved" expansion in powers of w. The exact value is .D(0.3) = 0.563683. The 
results close to the exact value are indicated with a greater number of digits. 
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N 


(a) 


(P) 


(c) 


(d) 


(e) 


(f) 


3 


-0.51 


1.04 


0.925 


1.01 


0.739 


0.97 


5 


-15 


1.9 


1.22 


1.05 


1.7 


0.8769 


10 


4 10 4 


-75 


2.4 


0.8412 


2.2 


0.857092 


12 


2 10 6 


349 


3.6 


0.678 


1.8 


0.853068 


15 


-9 10 8 


3 10 3 


8.2 


0.37 


1.6 


0.852614 


20 


8 10 


—1 10 


56 


—0.27 


2 10 


0.853263 


25 


-2 10 19 


1 10 8 


597 


-1.06 


0.59 


0.853463 


30 


1 10 25 


-3 10 13 


8 10 3 


-2.0 


1.5 


0.853438 



Table 2: The same as in Table 1, for a s = 0.5. The exact value is D(0.5) = 0.853427. 



N 


(a) 


(b) 


(c) 


(d) 


(e) 


(/) 


3 


-9.1 


1.9 


1.7 


1.34 


1.8 


1.00 


5 


-296 


6.7 


3.2 


1.13 


6.1 


0.61 


10 


8 10 6 


-381 


17 


-0.442 


13 


0.69528 


12 


9 10 8 


-2 10 3 


41 


-1.35 


-18 


0.6873 


15 


-2 10 12 


2 10 4 


178 


-3.0 


—22 


0.70106 


20 


1 10 19 


5 10 6 


3 10 3 


-6.2 


-7.1 


0.69741 


25 


-2 10 24 


3 10 12 


6 10 4 


-10 


21 


0.69580 


30 


4 10 30 


-3 10 18 


2 10 6 


-14 


-12 


0.696204 



Table 3: The same as in Table 1, for a s = 0.8. The exact value is .D(0.8) = 0.696408. 
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again in all three Tables, the columns (f), (d) and (b) have better convergence properties 
than the columns (e), (c) and (a), respectively. As concerns the point (iii), we see from 
Table 3 that the salutary effect of the optimal conformal mapping is most spectacular 
when the damping of the exponential function exp (—u/a s ) is the weakest, i.e., when a s 
has the highest value, a s = 0.8 in our case. Indeed, in this case, only the combined 
technique of the optimal conformal mapping and the explicit treatment of the branch 
points leads to numerical convergence, with steadily increasing accuracy up to A^ = 30 by 
the least. 

In the case of a stronger damping (Table 1), the role of the optimal conformal mapping 
combined with a careful regard to the branch point singularities is again important, but 
good results are obtained also by the other methods, the success varying with the per- 
turbation order A" used; see the different columns (b) - (f) of Table 1 at different values 
of A". The asymptotic superiority of the optimal mapping (columns (e) and (f)) emerges 
at very high values of N; this mapping supersedes the other methods and turns out to 
be the best at least from A^ = 20 on up to the highest value of A" shown in the Tables, 
N = 30. 

It is not excluded, on the other hand, that even the best series, column (f), will 
exhibit numerical indications of divergence at still higher orders; note that the singularities 
survive in some form because we (on purpose, in order to simulate real situations) had not 
completely removed them (see(ii)). Consequently, as the Borel integral path runs along 
the cut, which in the w plane is mapped onto the boundary circle of the convergence disk, 
no convergence is warranted even in the column (f). It was already pointed out that for 
lower values of a s , where the influence of infrared renormalons is more suppressed (see 
e.g. Table 1), results close to -D(0.3) are obtained even when other methods are used, see 
Table 1, sometimes only at lower values of N. 

To further illustrate the use of the optimal conformal mapping introduced above, 



we investigated it in a model proposed in adjusted to better simulate the physical 
situation. We assume the case when the Borel function is exactly given by 

~ D 2 
B true {u) = l + D 2 u + -V + [B{u) - £ b n u n ) . (44) 

^ n=0 

In this expression the parameters D 2 and D% are for the moment arbitrary and the numbers 
b n are the Taylor coefficients of the expansion of B around the origin: 



OO 



B(u) = J2b n u n . (45) 

n=0 

For B[u) we choose the expression 

s go(2) ,^ Mn+AMu J^ BM+jmu 

B <"> = (T^) + g (TTW + S ' (46) 

with 

8(-l)' +1 (3/ 2 + 6/ + 2) 



Mi) 



3 / 2 (/ + l)(/ + 2) 2 
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Mi) 



16(-l) m (/ + 3/2) 



3 / 2 (Z + l)Q + 2) 2 
B (2) = 

B Q (l) = -Ao(-l) ; I > 3 

^(Z) = -AH) ; Z>3. (47) 



The function B(u) coincides actually with the Borel transform in the large (3$ limit |T4| , 
with finite numbers Nyy and Nm of UV renormalons and IR renormalons, respectively. 
The meaning of the the model adopted above for B trU e is clear: it represents a function 



with the first 3 Taylor terms specified explicitly on the right hand side of (fPj), and the 
higher order terms arising from the sum of Nuv UV renormalons and Njr IR renormalons. 
The perturbative expression of this model is therefore 

B pert (u) = 1 + D 2 u + . (48) 

The expansion to the same order in terms of the optimal conformal variable w can be 
obtained easily using (j33|): 

B pert , w (u) = l + D 2 C lW + (D 2 C 2 + ^C 2 )w 2 . (49) 
We consider also, for comparison, the expansion in terms of the variable (|29"D used in |23[ : 

B peH , z {u) = 1 + D 2 C lZ + (D 2 C 2 + yC^ 2 , (50) 



with Cn defined in (p5|) . We introduce now the expressions B true , B pert , B peTt:W and B perttZ 
in the Laplace integral (|13|) and define the corresponding quantities D trU e , D pert , D perttW 
and D perttZ . The integrals defining D pert and D perttZ are well defined, and for D true and 
D P ert,w we adopt the princupal value prescription (^4|) . 
Folowing [^3| we consider the ratio 



Hw = -J) ZTT) ' (51) 

and the similar quantity H z . Clearly, the inequalities \H W \ < 1 (or \H Z \ < 1) are the 
conditions for the accelerated methods based on the conformal mappings w (or z) to be 



successful. As in [^3| we look for the domain in the plane D 2 , D$ for which the accelerated 
methods give better results than the usual perturbation theory. In Figs. 2 (a) and (b) we 
represent these domains, for the conformal mappings z and w, respectively (we used in 
this example the value a s = 0.27). 

As seen from Fig. 2, the domains are bands bounded by parallel straight lines, one of 
the lines passing through the origin. The slope is given by 



/ °° du e- u ^ [u-CtW- C 2 w 
f™due- u / a '[u 2 -Cfw 2 ] 
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(52) 



15 




-101234 -101234 



D 2 D 2 

Figure 2: Bands corresponding to \H\ < 1 for the conformal mapping z of |23|] (a) and 
the optimal conformal mapping w defined in the present paper (b). 



for the conformal mapping ( P7|) with the coefficients C n from (p4|) . A similar relation 
defines the slope in the case of the conformal mapping (|2"9"|). As seen from fl5*2] ) the 
slopes depend uniquely by the conformal mapping and not on the details of the model 
function. The numerical values obtained from flo"2| ) are s z = —0.60 and s w = 0.98. 

As concerns the intercept of the second line defining the allowed domain, it is given 

by 

T J™&ue-^\B{u)-Y.lb n u n } 
/ °°d M e-^[ M 2 -CV] ' 



and depends on the model function B(u). With our choice ([46]) the intercept was rather 
stable when increasing the number of terms Njjv an d Njr in the expansion. The results 
presented in Fig. 2 correspond to N uv = 4 and N IR = 6. 

Fig. 2(a) shows, as already remarked in |23|], that the conformal mapping ( p9f ) brings 
no improvement when the low order coefficients D 2 and D3 are both positive, as is the 
case of physical interest (see @). On the contrary, as shown in Fig. 2 (b), there are pairs 
of positive (D 2 , D 3 ) for which the optimal conformal mapping improves the ordinary 
perturbation expansion (the point of coordinates D 2 = 0.724, D 3 = 1.23 obtained from 
(^) and ( |TT| ) is actually close to the upper boundary of the domain in Fig. 1(b)). Therefore, 
even at very low orders (N = 3) an improvement can be obtained in principle by using 
the optimal variable. As the first coefficients have no alternate signs, this might mean 
that the first IR renormalon competes with the first UV one in contributing to these 
coefficients. 

We recall that in the last model we compared only the conformal mappings, without 
additional information about the nature of the first singularities of the Borel transform. 
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3.3 Determination of a s {m? T ) from r decay 

As a final application of the method we discuss the determination of the strong coupling 
constant a s {m 2 .) from the hadronic r decay width. It is known that the theoretical error 
is at present the dominant ambiguity in this determination, and the main source of this 
error arises from higher orders in perturbation theory. This makes the hadronic r decay 
a very suitable place to apply the technique of conformal mapping, which accelerates 
the convergence of the perturbative expansion and reduces the truncation error. As we 



mentioned, this problem was studied previously in p3[ , where the conformal mappings 
were used to reduce only the effect of the UV renormalons. It is of interest to use also 
the optimal conformal mapping, whose properties were demonstrated on mathematical 
models. We do not attempt to make here a complete analysis of a s (m^) determination, 
but only point out the effect of the combined technique of optimal conformal mapping 
and the dimplementation of the correct behaviour of the Borel transform near the first 
singularities. 

We used as starting point the Borel sum ( |15"D of R T and evaluated this expression using 
both the standard Taylor expansion (|ll|) of the Borel transform in powers of u, and the 
optimized expression (|30|) proposed by us. For comparison with previous work, we notice 
that the "standard expansion" in our approach is equivalent to the method of integration 



along the circle proposed in in the particular case of the one loop running coupling. 
The expansions were truncated at N = 2, with the coefficients h n determined from (|9]) 
and flT2"|). In the improved expansion ( |3"0"D we used the values 71 and 72 given in (^) and 



the coefficients c n were computed such as to reproduce the first three coefficients b n from 
(0)- 

In Fig. 3 we give the results corresponding to the standard Taylor expansion (|ll]) of 
the Borel transform (curve (a)) and the improved expansion ([It]) (curve (b)), for various 
values of a s (m^). Using the experimental value |31| 



{Rr)ex P = 3.645 ± 0.024 , (54) 

we obtain 

a s (m 2 T ) = 0.343 ±(0.009) exp , (55) 
using the standard Taylor expansion (|TT|) , and 

a 8 (m 2 T ) = 0.318 ± (0.007) e:Ep , (56) 



using the optimised expression (j30|). The improved expansion leads to a value of the 
coupling constant a s (m^) lower by about 8% than the result given by the standard Taylor 
expansion of the Borel transform. Actually, as in the above discussion of the model 
functions at low N and similar values of a s , the major contribution in shifting the value 
of a s towards smaller values is brought by the explicit treatment of the first singularities 
of the Borel transform. At the small values of a s relevant for the present problem, the 
effect of the conformal mapping is barely seen. 

In (|55| ) and (|56| ) we indicated only the experimental error, which is very small. On 



the other hand, it is not easy to ascribe a definite theoretical error to these results. The 
problem of the theoretical error of a s (m 2 ) was discussed in many papers, in particular in 
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3.9 



3.8 - 



3.7 - 




3.6 - 




3.5 



0.28 



0.3 



0.32 



0.34 



0.36 



0.38 



Figure 3: The Borel summation 
and the improved expansion ( |30[ 
experimental values. 



flU) 

(b), 



a s (m T ) 



of R T , using the standard Taylor expansion (a) 
as functions of a Ami). The band indicates the 



T7| , p3| , P5fl-P0[, with different conclusions about its magnitude. One can safely neglect 



the effect of the uncertainties in the QCD parameters (quark masses, gluon condensate 
etc), which is small [27], [[2JJ (leaving aside the still open problem of the 1/s corrections). 



The ambiguities related to the prescription chosen for computing the Laplace integral 
are believed to be small too, due to the conjecture that these ambiguities must be com- 
pensated by corresponding ambiguities in the condensates. The most important sources 
of theoretical error remain therefore those related to the analytic continuation from the 
euclidian to the minkowskian region, and the truncation of the perturbative expansion. A 
complete discussion of these errors is outside the objective of this paper. Concerning the 
analytic continuation, we only mention that in the derivation of ([15]) the perturbative ex- 
pansion of the Adler function was assumed to be equally valid in the euclidian region and 
in the complex plane near the timelike axis, which is certainly not true. It is not trivial 
to relax this assumption and see its impact on the determination of a s (m^). As concerns 
the truncation error, the estimate 8a s {m1) ~ 0.05 was suggested in and P5| , by com- 
paring the predictions of different summation procedures. In |Kj it was claimed on the 
other hand that much less errors are obtained if the renormalization group invariance of 
the perturbation series is exploited in an optimal way. The present work points towards a 



similar conclusion: indeed, as was remarked also in ]17[], it is rather arbitrary to interpret 
the spread of the results produced by different conformal mappings as a measure of the 
theoretical error, as suggested in [ffijfl . Our investigation on mathematical models shows 
that the truncation error depends on the choice of the conformal mapping, being smaller if 
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more information on the analyticity of the function is taken into account. The expansion 
proposed in our work exploits in an optimal way the (renormalization group invariant) 
information on the first renormalons of the Borel transform, and we therefore expect that 
the truncation error of the result (BH) is smaller than the estimate given above. 



4 Conclusions 

The technique of the optimal conformal mapping of the Borel plane, discussed in this 
paper, can be seen as an alternative resummation of higher-order effects in perturbative 
QCD. This resummation method has a physical content in the sense that the requirement 
of convergence in powers of the optimal variable w(u) amounts to a statement on ana- 
lyticity in the whole double-cut Borel plane. Indeed, the theorem [|22| on the asymptotic 



rate of convergence of power series, on which it is based, is dependent upon the condition 
that the function f(u) (which is expanded) and the function w(u) (in powers of which 
f(u) is expanded, see (pip) should have the same location of singularities. The method 
of the optimal conformal mapping allows us to make full use of this analyticity property. 

This remarkable feature is lost if the function is expanded in powers of some other 
variable, be it u or a conformal mapping of u such that only a part of the analyticity 
domain is mapped inside the convergence circle. In such cases, the convergence domain 
is smaller than the region of analyticity, and the requirement of convergence has to be 
supplemented with the analyticity condition. Only in the case of the optimal mapping 
the two regions are identical. 

As renormalons express the properties of the Feynman diagrams of the process, a 
statement about their location implies a statement about the physics of the process con- 
sidered. 

If the power expansion is truncated at a definite order, as is the case in practice, the 
roles of u and w(u) are modified. While a polynomial in u is holomorphic in u and has 
no singular ites in the Borel plane, a polynomial in w has the same analyticity region as 
the expanded function, having the cuts equally located. Since singularities have physical 
interpretation, every polynomial in w(u) carries this piece of information. 

We demonstrated the practical use of the optimized expansion numerically on a large 
number of model functions, for a sufficiently large truncation order (N > 5). The accuracy 
of the Laplace integral is especially increased by the optimal variable if the coupling 
constant is large and the exponential damping of the integrand is weak. In these cases the 
knowledge (even approximate) of the behaviour near the first renormalons, combined with 
the expansion in the optimal variable, leads to very accurate results, while the expansion 
in the Borel variable, though partially improved by the treatment of the branch points, 
fails dramatically. On the other hand, at low orders of perturbation expansion and for 
values of the coupling constant of physical interest the effect of the optimal conformal 
mapping is not very visible and the predominant effect is given by the explicit treatment 
of the nearest branch points. This was actually the case with the determination of the 
strong coupling constant a s (m*) from the hadronic r decay width: the combined technique 
of conformal mapping and the explicit treatment of the first branch points of the Borel 
transform reduce by about 8% the value given by the usual Taylor expansion in the Borel 
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variable. The major contribution to this result is brought by the theoretical information 
IIJ, about the nature of the first renormalons. 
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